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Abstract 

We build on existing calculations of the wall velocity of the expanding bub- 
bles of the broken symmetry phase in a first-order electroweak phase transition 
within the light stop scenario (LSS) of the MSSM. We carry out the analysis 
using the 2-loop thermal potential for values of the Higgs mass consistent with 
present experimental bounds. Our approach relies on describing the interac- 
tion between the bubble and the hot plasma by a single friction parameter, 
which we fix by matching to an existing 1-loop computation and extrapo- 
late it to our regime of interest. For a sufficiently strong phase transition (in 
which washout of the newly-created baryon asymmetry is prevented) we ob- 
tain values of the wall velocity, ~ 0.05, far below the speed of sound in the 
medium, and not very much deviating from the previous 1-loop calculation. 
We also find that the phase transition is about 10% stronger than suggested 
by simply evaluating the thermal potential at the critical temperature. 
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1 Introduction 



It has long been argued that Sakharov's three conditions for baryogenesis [1] (devi- 
ation from thermal equilibrium, CP violation, and baryon number violation) could 
have been met in the context of a first-order electroweak phase transition (EWPT) 
at the ~ 100 GeV scale [2] (for a review, see e.g. [3]). Baryon number violation in 
this setting would have proceeded through so-called 'sphaleron' transitions, which 
are Boltzmann-supprcsscd and become slow after the EWPT at T < 100 GeV, with 
the sphaleron energy Espi,/T v/T, v being the Higgs VEV. In a first-order transi- 
tion separate regions of the old and the new phase coexist, with bubbles of the new 
phase nucleating and expanding in the old phase until they fill all space. Baryon 
number violation would have taken place outside the growing bubbles, and then the 
baryon asymmetry would have been transported across the bubble wall into the new 
phase, where it must avoid annihilation by washout processes. This is granted by a 
sufficiently strong phase transition with v/T > 1.0 [4\, which is impossible to satisfy 
within the Standard Model for experimentally allowed values of the Higgs mass [5]. 
Therefore extensions of the Standard Model, such as models with extra Higgs fields 
[6,7] or singlets [8-13], or non-standard Higgs potentials [14,15] are required for 
successful electroweak baryogenesis. 

The possibility of a strong electroweak phase transition in a minimal super- 
symmetric setting leading to the production of a baryon asymmetry consistent with 
observations has been studied extensively. Electroweak baryogenesis has been shown 
to be feasible in a specific region in the supersymmctric mass parameter space. This 
setting is generally known as the light stop scenario (LSS), characterised by a (pre- 
dominantly right-handed) light stop with a mass lighter than or comparable to that 
of the top quark [16-19]. All other squarks and sicptons are typically taken to be 
at a much higher mass scale. Being one more light bosonic species (in addition 
to the weak gauge bosons) that couples to the Higgs, the light stop increases the 
upper bound on the Higgs mass compatible with a strong phase transition to about 
127 GeV [19]. In turn, the predominantly left-handed stop must be heavy to agree 
with electroweak precision tests and to provide a sufficiently heavy Higgs boson. 
Gluinos are generally considered heavy and thus decoupled from the thermal bath 
in order to suppress their potentially large contribution to the effective thermal light 
stop mass. Charginos and neutralinos should remain light as they provide the only 
additional CP- violating currents available in this context, needed to generate the 
observed baryon asymmetry of the Universe [20-23]. The CP-odd Higgs mass is 
large to avoid potentially large contributions to the electric dipole moments of the 
electron and the neutron, leaving one light, SM-like Higgs boson. 

A crucial parameter entering the computation of the generated baryon asymme- 
try [20-23] is the velocity of the expanding bubble walls, Vy^. In particular, it is 
essential that the wall velocity is smaller than the speed of sound in the plasma, so 
that diffusion of charges from the bubble wall into the symmetric phase is possible 
(see however ref . [24] for a recent proposal for supersonic electroweak baryogenesis) . 
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The wall velocity is determined by the friction induced by the motion of the 
bubble through the plasma, and by the pressure difference across the bubble wall. 
Microscopically, friction is related to deviations from equilibrium in the plasma. In 
a semiclassical approximation this effect can be described by a set of Boltzmann 
equations coupled to the equations of motion of the Higgs field. Such a system of 
equations was written down and solved for the first time in ref. [25] in the context of 
small Higgs masses in the minimal Standard Model. Later on the computation was 
generalized to the LSS in ref. [26]. While the typical wall velocity in the Standard 
Model was found to be Vy^ ~ 0.35, its value in the LSS was determined to be almost 
an order of magnitude smaller. This is due to additional friction induced by the 
light stops. 

Ref. [26] makes use of the thermal Higgs potential in 1-loop approximation, while 
2-loop corrections are known to be vital to obtain a strong phase transtion for realis- 
tic Higgs masses [16-19]. The aim of the present work is to extend the computation 
of the wall velocity in the LSS to the case of the 2-loop thermal potential, result- 
ing in a sufficiently strong phase transition for experimentally allowed Higgs masses 
around and above 115 GeV. 

Rather than repeating the microscopic analysis of friction from the literature, we 
model friction by a single friction parameter, which we fix by matching to ref. [26] and 
extrapolate it to our Higgs mass range of interest. As we will show such an approach 
leads to a much simpler set of equations, basically relativistic hydrodynamics coupled 
to the equation of motion of the Higgs field [27]. This method has been successfully 
used in ref. [28] to numerically simulate bubble growth. The main idea of our work 
is that friction is mainly related to interactions between the wall and the hot plasma, 
and does therefore hardly change when higher-order corrections are included in the 
thermal potential. On the other hand, friction would indeed change significantly, if 
the composition of the plasma would change, e.g. by removing the light stops. 

We will show that the main result of ref. [26] carries over to the 2-loop case: 
bubble walls in the LSS move very slowly, ~ 0.05, far below the speed of sound 
in the medium. This is also true for Higgs masses around and above 115 GeV, and 
for a phase transition strong enough to avoid baryon number washout. As a by- 
product we also arrive at a more reliable determination of the strength of the phase 
transition, reading off the Higgs VEV inside a realistic bubble, including effects of 
reheating, rather than taking it from the minimum of the thermal potential at the 
critical temperature. Thus we find that the relevant value of v/T is increased by 
about 10%. 
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2 Calculation of the wall velocity 



2.1 The perfect fluid setting and the hydrodynamic equa- 
tions 

We treat the problem of the moving bubble wall by modelling the hot plasma as 
a perfect relativistic fluid, assuming conservation of its energy-momentum tensor, 
which is the sum of the separate contributions from the fluid and the Higgs field 
[27]: 

+ {p + P)u''u'' - Pg""^ = 0. (1) 

Here p the energy density, P the pressure, and g^'' the usual Minkowski metric. 

Stationary^ solutions to the problem of the expanding bubble wall are usually 
divided into two categories depending on whether the wall advances at a velocity 
above or below the speed of sound in the medium [29-31]. In subsonic solutions 
('deflagrations') the bubble wall is preceded by a 'shock front' that accelerates and 
heats up the plasma, which is brought to rest by the bubble wall passing through. 
In supersonic solutions ('detonations') the plasma is hit by the bubble wall while 
at rest and accelerated, and is brought back to rest by a rarefaction wave which 
follows the wall. For the purposes of this work velocities remain well below the 
speed of sound and we will have no cause to explore supersonic solutions to the 
hydrodynamic equations. 

It is convenient to work with the stress-energy conservation equation in the rest 
frame of the expanding bubble wall, picking the radial direction and approximating 
the situation at the bubble wall as planar. We are left with one spatial coordi- 
nate and no time dependence. Applying the relevant thermodynamic relations and 
introducing our friction coefficient rj we arrive at the coupled system [27, 32] 

dM^) _ dV{^,T) ^ ^V\,^M^) .2) 
dx^ dip Tsi dx 

(4aT'-T?^)j\ = C. (3) 

(4«r* - T?^^^)j\' + - r) + = C2 (4) 

where (p is the scalar field, v is the fluid velocity, T the temperature, 7 the relativistic 
factor (1 — v^)~^/^, Tgi the plasma temperature in the symmetric phase ahead of 



^Steady-state solutions to the hydrodynamic equations for a spherical bubble are similarity 
solutions, that is, they maintain their relative shape but rescale as the bubble grows. 
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the advancing wall, and Ci, C2 are integration constants best determined in the 
symmetric phase where most contributions vanish. The radiative pressure is Pr = 
aT^ = ^g*T'^, g* being the number of effective degrees of freedom in the plasma at 
the temperature T. The special form of the friction term in (2) is motivated by the 
microscopic analysis of ref . [25] . 

This system can be solved numerically across the bubble wall by imposing van- 
ishing (f derivatives at both boundaries of the integration interval and vanishing 
Higgs VEV in the symmetric phase. We continue the computation of the bubble 
profile by solving the hydrodynamic equation in the region between the bubble wall 
and the shock front, for which we drop the planar approximation and take into 
account the sphericity of the bubble (note that on this step we are already in the 
symmetric phase where we take (/? = V{lp,T) = 0). Lastly we solve the conserva- 
tion equation across the shock front to relate our findings to the temperature of the 
universe unperturbed by the bubble [31]. 

As usual, a first-order phase transition is indicated by the effective potential at 
the relevant temperature, which has a global minimum at zero Higgs VEV, devel- 
oping a second local minimum at nonzero VEV as the temperature of the universe 
decreases (see e.g. [33]). As T keeps decreasing the value of V at the second minimum 
approaches the value of the symmetric minimum. The minima become degenerate 
at the so-called critical temperature, Tc. Nucleation of bubbles of the new phase 
becomes possible for T < Tc. The phase transition is deemed to begin once the 
integrated probability of bubble nucleation in the horizon volume 

P{T) = dP = {T/Vol)-VH-dt = jf^e-^'^^jr 

reaches unity, at the nucleation temperature T„ (here we have taken F/ Vol — 
A^(T)e~-'^'=/^ ^ T^e"^^/^, Fc being the free energy of the so-called critical bubble, 
just large enough to spontaneously grow without collapsing under surface tension, 
at the temperature T). Once bubbles begin to nucleate they grow and occupy the 
whole space very quickly, ending the phase transition. 



2.2 The wall velocity in the MSSM 

The wall velocity of expanding bubbles of the broken symmetry phase in a Ist-order 
electroweak phase transition was calculated microscopically for the Standard Model 
by Moore and Prokopec [25] and for the MSSM, based on the same procedure, by 
John and Schmidt [26]. The equation of motion for one scalar background field 
coupled to the distribution functions of the particles in the plasma is [25] 

dV((fi) ^dm'^ f d^p . ^ 

Here V{(f) is the renormalised vacuum potential, the sum is over all particle species, 
and the mass dependence of each particle on the Higgs VEV is given by its respective 
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couplings. The vacuum potential combines with the equilibrium part of the particle 
distribution to produce the finite-temperature effective potential Veff{(p,T). The 
equation of motion can then be written as 

where 5f expresses the deviation from equilibrium of the relevant particle popu- 
lations, which is responsible for the wall friction. The evolution of the particle 
populations, in the semiclassical approximation, is given by Boltzmann equations of 
the form 

^^=dtf + ^-dgf+p- dpf = -C[/]; (8) 

C[f] being the collision integral. These Boltzmann equations together with the 
equation of motion for the scalar field constitute a very complicated set of coupled 
integro-differential equations, which is very difficult to solve. 

Several simplifications are possible in order to solve the equation of motion. 
Taking C[f] = constitutes the 'free particle' approximation. In this approach the 
Boltzmann equation can be solved exactly [25] and / integrated over momentum. 
Recently this approximation was used in ref. [34] to show that under certain circum- 
stances the bubble wall velocity may approach the speed of light. In the 'relaxation 
time approximation' C[f] = (— ), with r independent of momentum. This approach 
becomes particularly simple if t » L^j with the wall thickness, because then 
the spatial derivatives of 6 can be neglected. 

The approach first followed by Moore and Prokopec for the Standard Model [25] , 
then by John and Schmidt for the MSSM [26] (the 'fluid approximation') describes 
/ for each relevant particle as an equilibrium distribution where E/T is modified by 
perturbations, 

where 5 — —[SjA, + 6jj,bg + ^{5T -\- ST^g) + Pz{5v + 5vi,g)]. It is sufficient to treat 
in this way only 'heavy' species which couple strongly to the Higgs (top quarks, 
W bosons, and Z bosons for the Standard Model, and additionally right-handed 
stops in the MSSM). Note that in the MSSM, as mentioned, other ('superheavy') 
particles are taken as decoupled from the plasma and show up only in the shape 
of the effective vacuum potential of the theory. The treatment of 'hght' particles 
is simplified by treating them as one species with common perturbations to the 
equilibrium chemical potential, temperature and fluid velocity S^bg-, SThg and Svjjg 
(as opposed to the individual perturbations for 'heavy' species 6fi. ST and Sv). 

With this approximation the Boltzmann equation for each particle species can be 
expanded to linear order in the perturbations and then J d^p/{2n)^, J Ed^p/{2TT)^, 
and J Pzd^p/ (27r)^ integrations carried out. This, plus the fluid equation, constitutes 
the system to be solved numerically ^. 



*John and Schmidt, in the MSSM case, refer expUcitly to the simphfication in which all S' = 0. 
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3 The MSSM 2-loop potential 



As the model's effective tlieory we base ourselves on the 2-loop, finite temperature 
MSSM potential calculated in refs. [16-19, 35] with only one (hght) background field 
with SM-like couplings to vector bosons and fermions. Only third generation squarks 
are considered to be at the clectroweak scale, and we assume no mixing between 
right- and left-handed stops. With this the additional supersymmetric parameters 
of the model are just the soft supersymmetry-breaking mass parameters for the 
left- and right-handed stop mg and niu and tan/3 = ^ (with vi, V2 being the zero 
temperature expectation values of the real parts 0i, 02 of the neutral components 
of the two supersymmetric Higgs doublets). In the present case of only one light 
Higgs field the effective potential can be expressed, as in the Standard Model, as a 
function of only one background field. 

3.1 The 1-loop potential 

The 1-loop portion of the effective potential for the light Higgs field ip is [35] 

Vtreei^) = -^/^V + ^fV COS^ 2/3 (10) 
T4-1oop,0-t(^) = E ^2^tm^S ^ - a] (11) 
l^l-loop,thermal(¥^, ^) = E ^^^4^^^^] (12) 

Note that we take g' = 0. Since 1-loop contributions to the potential are comparable 
to the tree-level portion, the parameter //^ in the tree-level part is chosen so that 
the minimum of the total 1-loop, non-thermal potential lies at cpo — 245.7 GeV. 

Sums run over all species that contribute significantly. For the 1-loop part this 
includes stops, tops and W and Z bosons. The left-handed stops do not contribute 
to the thermal piece. The number of degrees of freedom for each species, rii, is 

nt = -12, = = 6, n^y = 6, = 3 (13) 

We take Q = mz and Cj = | for vector bosons, | for scalars and fermions. The 
relevant expansion of the J functions in the 1-loop, thermal bit is of the form 



bosons \jk)= dx X log \1 - exp + (^) 




In that case a fluid equation for eacli scalar field formally identical to (2) is produced with a 
calculable friction parameter r]. 
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J fermior 

= -(27r^) 

The notation m,j (for bosons) in the 1-loop, thermal piece indicates resummed 
masses, m,^ are obtained from by adding the leading temperature-dependent 
self-energy contributions (see below) . An effect of the 2- loop calculation is to resum 
all masses in the thermal piece of the 1-loop potential^. 
The expressions for the particle masses are: 

<pM = ^^tV'sm'/3, (16) 

^iM) = ^1 + + dU<p), (17) 

mU^) ^ml + m'M + D}^{<p) (18) 
assuming no mixing between left- and right-handed stops and taking 

= ''''' ^^^^ 

DU^) = \{^ sin^ Ow^ ^ V cos 2/3 (20) 

Note that with the choice g'' = we have = for all </?. 

As field-dependent mass for the W and Z gause bosons we take 

m^((^) = m|((^) = ^^V^- (21) 
3.2 The 2-loop contributions 

Following the notation of ref. [35] , the relevant Standard Model 2-loop contribution 
to the potential can be written as 



T/(2) _ rj.2 



\.-2f^, Ml 51, M\ 3. ,2 .,.2N, M + 2Ml " 
^ [l ^ " ¥ ^ Y j + 2 (^"^ " ^""^ 37^ + ^^^^ 



647r2 



16ff2 { -log2---CB 1 + 9ht sin' /3 I -log2- CB 



(22) 



^In the 1-loop version of our potential (which we use to reproduce the results of ref. [26]) 
the masses in the 1-loop, T-dependent piece are not resummed. Following ref. [35] we do re- 
sum the bosonic masses (only the longitudinal degrees of freedom for the gauge bosons, pho- 
tons included) in the term cubic in m in the expansion, through the addition of the piece 
AV{(fi,T) = —^:^T,ni[fhf{(l),T) — mf{ip,T)] running over the relevant species with uwl = 2, 
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where cb — log(47r) — 7^, 7^ 0.577215665 being Euler's constant. We take the 
strong couphng and top Yukawa couphng respectively as gs ~ 1.228 and ht = . 
Here Af^ = \g^(f^ is the weak gauge boson mass and M| = + '^g^T'^ the 
longitudinal resummed mass (corrected by the self-energy). 

We have retained the supersymmetric contributions which relate to non-decoupled 
species in the plasma. The relevant diagrams are shown in [35] . We take as the total 
supersymmetric, 2-loop part of the potential for our calculation 



V2-I00P, MSSM — 



167r2 



167r2 



^in log 



3T 



+ 



32.^ (A^-^/J) log 3^ 



*'^'(iV=-l)fe-l)m? + 



tR 



647r2' " 



9j 
6 



{N^ + l)mf + -h'^t sin^ /3 {mhrrii^ + Sm^mt^) 



(23) 



We take the number of colours as 3, and C2 ~ 3.3025. Note that, given the 
large 1-loop corrections to the Higgs mass, rfih and are expressed on the basis 

on taking ruh = — ^"2'°°'' at the minimum of the 1-loop potential. 

The region in rriij-mQ space for tan (3 — 4 which provides both a sufficiently 
strong phase transition and presently acceptable values of the Higgs mass is shown 
in figure 1. We observe that a strong phase transition, i.e. ,^ = |f > 1, can be 
achieved for a light stop mass around 135 GeV and have stop mass above 10 TeV. 
We will see later on that this procedure underestimates the strength of the phase 
transition by about 10%. Here the Higgs mass is around 115 GeV. Somewhat larger 
Higgs masses can be obtained by using larger values of tan (3 and/or rriQ. 



4 Calibration of the friction parameter 

Our goal is to calculate the wall velocity in the region of parameter space which 
provides a sufficiently strong phase transition and is compatible with present bounds 
on the mass of a SM-like Higgs boson. We assume that the phenomenological friction 
parameter 77 depends exclusively on mu and tan (3, as these are the parameters which 
govern the abundance and coupling of the light stops. We match the 1-loop results 
of John and Schmidt through the system (2)-(4) to determine our friction parameter 
rj for tan (3 e [2, 6], mu G [—60, -f-60] GeV. The coupled system can be linearised and 
solved numerically. The fitted mu, tan /3- dependent friction parameter can then be 
extrapolated to the relevant region in MSSM parameter space. The values of the 
fitted friction parameter which reproduce the results of John and Schmidt are shown 
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Figure 1: Values of the strength of the phase transition ^ = |f and the Higgs mass 
in the region of parameter space (with tan ^ = 4) of interest for baryogenesis. For 
this value of tan ji the right-handed stop mass (we assume no mixing) varies in this 
range from m^^ ~ 135.2 to m^^ ~ 138.8 GeV and the left-handed-stop mass is given 
by mq in units if GeV. 
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Table 1: Phenomenological friction coefficient 7] fitted to the 1-loop MSSM wall 
velocity in ref . [26] for mq = 2000 GeV and different values of tan /3, mu [GcV] . 



tan /3 


mfj 


(John and Schmidt) 


Fitted 


2 


-60^ 


0.060 


4.58 







0.090 


3.36 




+60^ 


0.160 


1.92 


4 


-60^ 


0.080 


4.35 







0.115 


3.16 




+602 


0.140 


2.72 


6 


-602 


0.085 


4.65 







0.120 


3.06 




+602 


0.155 


2.55 



in table 1. We find that r] is almost independent of tan^, but grows significantly 
with lower stop masses, i.e. lower m^. The linear extrapolation to the mu interval 
of interest for tan/3 = 4,6 is shown in figure 2. As shown, the extrapolation to 
the required lower values of introduces some uncertainty which we quantify 
later on. Comparing with our treatment of SM-like friction in ref. [32] we find the 
friction coefficient enhanced by a factor 5 to 10. This seems very plausible, given 
the couplings and number of degrees of freedom of the light stops compared to e.g. 
that of the W-bosons. 

5 Results 

With the extrapolated values of the friction parameter the coupled hydrodynamic 
equations (2)-(4) can be solved using the full 2-loop thermal potential of the LSS. 
The results for the wall velocity of this calculation are shown in table 2 in two rep- 
resentative cases for tan /3 = 4, 6 (these values being the most favourable choices as 
regards obtaining both a conveniently high Higgs mass and a strong phase transi- 
tion) . 

To test our extrapolation of the phenomenological friction parameter we found 
the highest possible rj variation found by taking only two calibration points from 
[26] for the case tan 13 = A. We show this two-point calibration in figure 2 and give 
the alternative wall velocities thus found for the case rrih = 112 GcV, fc = = l, 
in table 2. As seen, the wall velocity varies accordingly from 0.037 to 0.057. While 
this change is noticeable, the prediction of a very subsonic wall velocity in the LSS 
is robust, and justifies the use of ~ 0.05 in computations of the resulting baryon 
asymmetry, e.g. in ref. [22]. Keeping our original (3-point) calibration and sampling 
the most promising areas of parameter space as regards baryogenesis and the Higgs 
mass we find that the wall velocity varies less than 10% across the region shown in 
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-12000 -10000 -8000 -6000 -4000 -2000 2000 4000 

Figure 2: Linear extrapolation of 1] values to the mfj region of interest for tan /3 = 
4 (solid line, calibration values from [26] as triangles), tan/3 = 6 (dashed line, 
calibration values as crosses), m^^ G [135.2,184.1] GeV (almost no dependence on 

tan/3) 
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Table 2: Results of wall velocity calculation at 2-loop order for Higgs mass rrih = 112 
GeV, tan^ = 4,6 and strength of the phase transition ^ close to 1. The values of 
^f,— ^ in the broken symmetry phase are given for comparison. 



tan /3 


d — fc 
?c - 




Fitted f] 




4 


0.9 


1.01 


5.94 


0.044 




1.0 


1.14 


6.05 


0.043 






1.14 


4.58 (2-point) 


0.057 






1.14 


7.02 (2-point) 


0.037 


6 


0.9 


1.00 


6.66 


0.039 




1.0 


1.14 


6.79 


0.038 



figure 1, staying close to = 0.05. 

Relevant for baryon number washout is the value of the Higgs field inside the 
bubbles, Vb, related to the temperature inside the bubbles, T^. The resulting ratio 

— Vb/Tf, will be different from the commonly used = Vc/Tc, obtained from 
the equilibrium thermal potential at the critical temperature. It is important to 
understand that the determination of requires the solution of the bubble evolution 
equations in the presence of a plasma, such as eqs. (2)-(4). We show ^ in figure 3 
as a function of for the case tan^ — 4, mg = 14,000. is related to T^, 
while is related to Tf.. Note that in this case the nucleation temperature, T„, 
and the temperature inside the bubbles happen to be nearly identical. The value 
of is significantly higher than ^c, so the commonly used criterion systematically 
underestimates the strength of the phase transition by at least 10%. 

The low values of the wall velocity found in this setting question some of the 
assumptions often made when studying the real time history of first-order phase 
transitions. The bubbles' relatively slow expansion increases the gap between the 
nucleation temperature and the so-called finalisation temperature Tf, at which the 
proportion of space occupied by the growing bubbles reaches unity and the phase 
transition ends. As an approximation it is usually assumed in the calculation of Tf 
(see, e.g. [33]) that the bubbles expand at the speed of light. It seems reasonable in 
our case to quantify the error resulting from that assumption. An alternative choice 
of e.g. Vw = 0.05 results in a different finalisation temperature Tj. As an extension 
to our calculation we investigate the variation in the wall velocity introduced by 
assuming that the temperature of the undisturbed universe outside the bubble is 
a) the nucleation temperature T„; b) the finalisation temperature calculated in the 
usual way Tf] and c) the finalisation temperature calculated assuming a slow wall 
velocity Tj.. We show the results in table 2 for the case mh = 116 GcV, mfj = —11900 
GeV^ 6 with tan /3 = 4, 6. The choice of finalisation temperature does not introduce 
a huge variation in but whether we take the temperature of the undisturbed 

^Low values of mfj are most favourable for baryogenesis. 
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0.85 
-12000 



-11500 



(GeV^) 



Figure 3: Values of = v/T for the case tan/3 = 4, mq = 14, 000 GeV {rrih ~ 112 
GeV) at the critical and nucleation temperatures and in the broken symmetry phase 
where sphalerons must be suppressed to avoid washing out the newly-generated 
baryon asymmetry 



13 



Table 3: Wall velocity for the case rrih = 116 GeV, = —11900 GeV^ taking as 
the temperature of the universe the nucleation temperature and the finalisation 
temperatures calculated on the assumption Vyj — 1 {Tf), — 0.05 (T^) 



tan/3 T i = Fitted 77 Vw 

~ % 0^92 6l0 0.044 

Tf 0.94 0.051 

r; 0.94 0.054 

6 T„ 0.95 6.88 0.039 

Tf 0.96 0.045 

T'f 0.96 0.047 



universe as T„ or Tf makes a significant difference. So the bubbles accelerate between 
the time when they are nucleated and the time when they collide at the end of the 
phase transition, which is the relevant time for electroweak baryogenesis. However, 
it is important to recall that our hydrodynamic treatment assumes bubble expansion 
at a steady-state velocity and cannot account for bubble wall acceleration, due to the 
variation in temperature between the beginning and the end of the phase transition 
or for any other reason. 



6 Conclusions 

We have computed the bubble wall velocity during a strong first-order electroweak 
phase transition in the light stop scenario of the MSSM. This requires a Higgs mass 
close to the present lower experimental bound of 115 GeV, a light right-handed 
stop with a mass around 140 GeV and a heavy left-handed stop with a mass in the 
multi-TeV range. 

The bubble wall velocity is determined by balancing the pressure inside the 
expanding bubbles with friction from the hot plasma. We describe the system as 
a single Higgs field coupled to relativistic hydrodynamics via a phenomenological 
friction term. The size of this friction term we determine by matching to the wall 
velocities computed microscopically in ref. [26]. Ref. [26] used the thermal Higgs 
potential in the 1-loop approximation and was therefore not able to have a strong 
phase transition at realistic Higgs masses. We use the full 2-loop potential and are 
able to treat the physically interesting case oi v/T ^ 1 for Higgs masses of around 
and above 115 GcV. We find a wall velocity ~ 0.05, not very different from the 
results obtained in ref. [26] for much weaker phase transitions. The obtained friction 
coefficient is roughly 5 to 10 times larger than what is found in an SM-like situation 
[32] . It could also be used to treat e.g. the case of an additional singlet scalar field 
being added to the light stop scenario. Our result confirms the expectation of slowly 
moving bubble walls in the light stop scenario, in contrast to ref. [36], where Vy, ^ 0.4 
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was found, using a simplified microscopic model of friction. 

Having computed the full scalar profile of the expanding bubble, wc can more 
reliably check the criterion for avoiding baryon number washout, by taking the 
true Higgs field value and temperature inside the expanding bubble, rather than 
relying on the equilibrium potential at the critical temperature. This way the phase 
transitions turns out to be about 10% stronger. 

The wall velocity is not only crucial for the generated baryon asymmetry, but 
also for the production of gravitational waves (see, e.g. [37-39]). However, the 
gravitational wave signal goes down with the wall velocity, and should therefore not 
be observable for the MSSM electroweak phase transition. 

Finally, we will treat the non-supersymmetric case in more detail in ref. [40]. 
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